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1. Introduction 

The nonequilibrium work theorem is an equation in statistical mechan- 
ics that relates the free energy difference AF to the work W carried 
out on a system during a nonequilibrium transformation. The identity 
appeared in different, but as we show below equivalent, formulations 
in [1] and [2] in 1997, and in the series of papers [4]-[7] in 1977-1981. 
In physics literature, the identity is usually written in the form 

{e-P^) = e-/^^^, (1) 

where the average is taken over all possible system trajectories in 
the phase space, and /? is an inverse temperature. The identity first 
appeared in this form in [1] and [2]. Traditional equilibrium thermo- 
dynamics tells us that {W) ^ AF while the transformation of the 
system is infinitely slow. The identity (1) is a stronger statement, and in 
addition to this, it is valid for arbitrary transformations of the system. 
The identity is used effectively in computer simulations, as well as in 
experimental physics, to calculate the free energy difference between 
two states of the system by running many trajectories and taking the 
average value of e~^^ (see [8]~[12], [14]- [16] and references therein). 

The paper [3] discusses the connection between two different ver- 
sions of the identity, and shows that the papers [4]-[7] use a different 
definition of work. The identity obtained in [4]-[7] (referred to below 
as Bochkov-Kuzovlev's identity) reads: 
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where Wq is the work (in Bochkov-Kuzovlev's sense) performed on 
the system, and the angle brackets have the same meaning as in (1). 
The present paper shows that Jarzynski's and Bochkov-Kuzovlev's 
identities easily follow from each other. 

Since the identities involve taking an "average over trajectories", 
it is natural to interpret this average as the expectation relative to a 
probability measure on trajectories, while assuming that the system 
evolves stochastically. In terms of expectations the identities can be 
represented by the formulas 

£[6-^^^] = e-^^^ and E[e-^^°] = 1, 

where E is the expectation relative to a probability measure on phase 
space paths. For this probability measure, some analytical assumptions 
under which the identities hold are found. 



2. Notation and assumptions 

Let us assume that the evolution of our system is described by a Markov 
process Tt^Lo), t G [0, T], given through its transition density function. 
Let X = M^*^ be the phase space for our system, i.e. the set of values of 
r^. We assume that at time t = the distribution on the phase space 
X is given by the following density function: 

where ■ , A) : X — > R is a Hamiltonian parametrized by an externally 
controlled parameter A G A, A C is an open set; f3 = 1/(/cbT), 
fcs is the Boltzmann constant, T is the temperature of the system, 
-^Ao = Ix e~l^'^^^'^°^ dx is the partition function. We assume that for 
all Ao G A, 6"^^*^^''*'°^ dx < oo. We consider the situation when the 
external parameter A is a function of time [0, T] A, i.e. we actually 
consider a time-dependent Hamiltonian H{x,X{t)). Let E;^^ denote the 
expectation relative to the measure qxQ{x)dx. Below we assume that 
the changing in time external parameter A belongs to the space 

V[0, T] = {X: [0, T] ^ A, A = Ac + A,tep; Ac G C^, A^tep e C,tep} 

where = C^[0, T] is the space continuous function of bounded 
variation on [0,T], and 

Cstep = >Cstep[0, T] = {A( • ) = E-Jn ( " ) + An} 
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is the space of right continuous step functions corresponding to different 
partitions V = {0 = to < ■ ■ ■ < tn = T} of [0,T] and different finite 
sets of values {Aj}. 

To emphasize the fact that the function A G V[0, T] gives rise to the 
process Fj, we will use the notation T^. Let px{s,x,t,y), s,t e [0, T], 
s < t, X, y ^ X, be the transition density function for T^. Let X^'^''^^ 
denote the space of all paths [0, T] — > X. In terms of px we can construct 
a probability measure L;^ on X^^'"^ by means of the finite dimensional 
distributions 



f{u}{tQ),Lu{ti),. . . ,uj{tn))'Lx{dLo) = dxoqx{to)ixQ) (2) 

jdxiPx{to,Xo,ti,Xi) ■ ■ ■ jdXnP\{tn-l,Xn-l,tn,Xn) f{xo,Xi, . . . 

X X 

where {0 = to < ■ ■ ■ < 1,,, = T} is a partition of the interval [0, T], and 
/ : X"^^^ — > M is a bounded and measurable function. By Kolmogorov's 
extension theorem, the right hand side of this equality defines a prob- 
ability measure on the minimal c- algebra of X^'^''^^ generated by all 
cylindrical sets. We denote this a-algebra by adX^^''^^). We assume 
that px{s,x,t,y) satisfies one of the assumptions (see [17], chapter 2, 
paragraph 1) that guaranties that the measure L;^ is concentrated on 
the right continuous trajectories without discontinuities of the second 
kind. Also, we assume that the a-algebra cTc(X[0'^1) is augmented with 
all subsets of L^-null sets. 

We take X^^''^'^ as the probability space, i.e. we set f^ = X[0'^]. Then, 
L;^ is the distribution of F^, and L;^-a.s., F^(a;) = 

For each fixed A G A, we introduce another transition density 
function p{s,x,t,y,X) which represents the situation when the sys- 
tem evolves being controlled by a constant in time parameter. The 
following below Assumption 1 is the key assumption under which the 
nonequilibrium work theorem holds. 

Assumption 1. If X\^s,t] = then px{s,x,t,y) = p{s,x,t,y,\), where 
p{s, x,t,y,X) conserves the canonical distribution on the phase space 
X. Specifically, it satisfies the identity 



L 



Qxix) p{s, x, t, y, A) dx = qx{y). (3) 

X 

We make two further assumptions: 

Assumption 2. //A is constant on [s,t), and discontinuous at the 
point t, then for Lebesgue almost all x E X, 

lim / dyp{s,x,t- 6,y,X{s)) dz px{t - 6,y,t, z)f{z) 
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= dyp{s, X, t, y, A(s))/(y), 

where / : X — > R is bounded and continuous. 

Assumption 3. For all t € [0, T), and for all compacts K c X, 

lim sup / p{t,x,t + 6,y,X{t))f{y)dy - f{x) 
where f : X is bounded and continuous. 



0, 



Lemma 1. Let V = {0 = to < ti < ■ ■ ■ < tn = T} be a partition 
of [0,T], and let A( • ) = E"Jo^ ^i\ti,ti+i)( ' ) + ^n- We assume that if 
X\[s,t] = ^7 where X E A is constant, then px(s,x,t,y) = p(s,x,t,y,X). 
Further we assume that the function p\ satisfies Assumption 2. Then, 
for all s, t, ^ s < t ^ T , for all x & X, and for all bounded and 
continuous functions f : X ^ W, 

/ dypxis,x,t,y) f{y) = dxip{s,x,Ti,xi, X{s)) 

JX JX 

/ dx2p{Ti,xi,T2,X2,X{Ti))--- / dy p{Tk, Xk,t, y , X{Tk)) f {y) , (4) 

JX JX 

where {s < ti < ■ ■ ■ < Tk < t} = {V[j{s, t}) n [s, t] . Moreover, the right 
hand side of (4) defines a transition density function. 

Proof. First we prove that the right hand side of (4) defines a tran- 
sition density function. We have to verify the Chapman-Kolmogorov 
equation. 

Let r G (s, t), and V = VU{s,r, t}. Further let Vi = {s < n < T2 < 
■■■<Ti} = Vn[s,r], andT'a = {r < a < 6 < • • • < U < = rn[r,t]. 
Then, V1UV2 = V, and [r;, ^1) is the interval of the partition V such 
that r G [ti,^i). We have: 

/ dxlp{s,X,Tl,Xl,X{s)) / dX2p{Tl,Xl,T2,X2, X{ti)) ■ ■ ■ 
JX JX 

/ dyp{Ti,xi,r,y,X{Ti)) / dzip{r,y,^i,Zi,X{r)) ■ ■ ■ (5) 

JX JX 

I dZffi p(^m— 1) •2^m— 1 ) Cm) ^^mj -^(Cm— l)) P{Cmj •^mj Z, A(Cjn)) • 
JX 

Note that A(r) = X{ti), and hence, by the Chapman-Kolmogorov 
equation for the transition density function p( ■ , ■ , ■ , ■ , A(r;)), we have: 



Ix 



dyp{Ti,xi,r, y, X{Ti))p{r, y, ^i, zi, A(r)) = p{ti,xi,^i, zi, X{ti)). 
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Let us take a bounded and continuous function /, and show (4). 

/ „(.,x,*,,)/fe)d, = /dx;p(,,.,n-i.,x;,A(»)) 

Jx Jx 

/ dxipxin - Si,x'i,Ti,xi) 
■'^ (6) 
dx'kP{rk, Xk, t - 6k, x'k, A(rfc)) 



X 



dypx{t - Sk,x'k,t,y)f{y) 

which holds for all smaller than the mesh of V. Taking 

the repeated limit lim^^^o li™<52^0 " ' ' li™<5fc^0 of the right hand side 
of (6), applying Lebesgue's theorem, and taking into consideration 
Assumption 2, we obtain the right hand side of (4). □ 



3. Jarzynski's identity 

Let A = Ac + Astep be the decomposition of A G V[0, T] into a sum of a 
Ac G C^[0,T] and a Astcp G £stcp[0,T]. Further let dxTi : X x A ^ 
denote the partial derivative with respect to the second argument (i.e. 
with respect to the control parameter) . We assume that the restriction 
of dxTi. to X X A([0, T]) is bounded and measurable. Also, we assume 
that for each fixed A € A([0, T]), 7Y(-,A) is bounded and measurable. 
Everywhere below, the probability space O is the space We 
define the work Wx : — > R performed on the system by the formula 

Wxiu) = r{dxniT^iuj),Xit)),dXS))Ri 

^'n (7) 
+ ^(7^(r,^^(a;),A(t,))-W(r,^^H,A(i,-0))), 
1=1 

where {0 = to < ti < ■ ■ ■ < tn = T} are discontinuity points of A, ( , 
is the scalar product in D A, and the integral on the right hand side 
is the Lebesgue-Stieltjes integral, i.e. the sum of the Lebesgue-Stieltjes 
integrals with respect to the components of Ac. In the following, we skip 
the sign of the scalar product in the first term of (7), and simply write 

Wxiuj) = r dxn{r^{u;),xit))dx,{t) 

Jo 

n 

+ J2inirl{iv),xiu))-n{T^^{iv),xiu-o))). 

i=l 
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Let Fx = InZx (free energy of the system), and let AF = Fx(t) — 
-^A(o) (frse energy difference). Let El^ denote the expectation relative 
to the measure L;^. 

3.1. Jarzynski's Identity, case A G £step[0,T] 
Clearly, if A = Er=o ^''^lU,U+i) + ^ £step[0,r], 

n 

1=1 

Theorem 1 (Jarzynski's identity: case A G >Cstep[0, T]). Let A G 
Cstepi'^jT], and let the transition density function px satisfy Assump- 
tions 1 and 2. Further let ?{{ • , A) be bounded and measurable on X for 
each fixed A G A([0, T]). Then the function e~^^^ is "Lx-integrable, and 

Proof. L^-a.s., 

n 

Wx{o:) = Y.(H{u{ti), \) - H{uj{ti), Ai_i)) . 

i=l 

Without loss of generality, we can assume that t^ = 0, tn = T, adding 
these points with zero jumps if necessary. Note that Wx{uj) is a cylinder 
function. By Lemma 1, 

EUe-f'^^] = [ e-'3Sr=i(«(-(*«)'^')-^(-(^^)'^-i))LA(du;) 

= / dxoqxoixo) / dxip{to,xo,ti,xi,Xo) 
Jx Jx 

■■■ / dxn-ip{t n—2i Xn—2i tn—li ^n—lj 

Jx 
Jx 

(8) 

Note that Assumption 1 implies that for all s < for all y G X, and 
AG A, 

^^^ e-/3>^(,,A) Pis,^,t,y,X) = l. (9) 
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Taking into account this, and changing the order of integration in (8), 
we obtain 



Z\o Jx 

p g— /3?i(x„_i,A,i_i) 

dXn-1 ^-i3n{xn,\n-i) Pi*n-l,Xn-l,tn, Xji, K-l) ^^^^ 



^^'^ -pn(xiM) P(*o,2:o,ii,2:i, Ao). 



Starting from the end, we replace each integral in (10) with 1, which 
is valid by the relation (9), until we reach the very first integral (taken 
with respect to x„), which we replace with Z\^. Noticing that Aq = A(0) 
and An = A(r), we obtain 



A(0) 



z 



The theorem is proved. □ 

3.2. Jarzynski's Identity, case: A G C^[0, T] 

Theorem 2 (Jarzynski's identity, case: A G C^[0,T]). Let A G 

C^[0,T], and let the transition density function p\ satisfy Assumptions 
1, 2, and 3, and the probability distribution La of be given by (2). 
In addition, let the following assumptions be fulfilled: 

4- If G Cstep[0,T], and as n ^ oo, A" — >■ A uniformly on [0,T], 
then L^n — > weakly relative to the family of bounded continuous 
cylinder functions; 

5. The function 7{{-,X) is bounded and measurable on X for each 
fixed Xe X{[0,T]); 

6. The function d\H is bounded on X x A([0, T]); 

7. The functions d\H{x, ■ ) are equicontinuous as a family of functions 
parameterized by x E X; 

8. The function d\H{ • , A) is continuous for each fixed A G A. 
Then, the function e~^^^ is 'hx-integrable, and 

Ei^Je-/^^-] = e-/^^^. (11) 
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Lemma 2. Let X e C^[0,T]. Then, under Assumptions 6-8, the 
function W\{u) is Uc-measurable. 

Proof. Below, for an arbitrary small e we construct a a"c-measurable 
function F : O — > M such that 

snY>\F{io) -Wx{io)\ <e. (12) 
n 

We find a function Agtep = TJ^Zq KU)\ti,u+-L) + Ktn) e >Cstep[0,r] such 
that for all a; G for all t G [0,T], 

\dx'H{T^{u:),\{t)) dxn{T^{u:),\,,^{t))\ < e. 

This is possible by Assumption 7. Thus, it suffices to find a measurable 
function F verifying (12) for Wx of the form fldxn{T^{uj),\)dXc{t), 
where A G A([0,T]) is fixed, and Ac is not a constant on [s,r]. 

By assumption, L;^-a.s., the paths of are right continuous and 
without discontinuities of the second kind. This implies that the map 
: [0,r] X O ^ X is (^([0,r]) O (7c,'^(^))-measurable (where 
^([0,T]) and 3S{X) are Borel a-algebras) which follows from [17] 
(Chapter 2, Theorem 11). By Assumption 6, d\7i{ • , • ) is bounded, say 
by a constant M. We divide the ball of radius M in M' with the center in 
the origin, by a finite number of sets Oj whose diameter is smaller than 
e/Vg^[Ac], where Vq^[Ac] denotes the variation of Ac on [0, T]. Further 
let Ai = dx'H{-,X)-\Oi), and Q = {r^)-\Ai) c [0,r] x fi. The 
sets Ai are open by Assumption 8. The sets Cj are ^([0, T]) ctc- 
measurablc by the (^([0, T])(E>crc, =^(X))-mcasurability of the map F'^. 
Fix Xi G Ai, and consider the function ^{uj, t) = J2i 9x'H{xi, A)!^. (w, t). 
Clearly, suptg[o,T],c^€0 |5AW('^(i), A) - $(u;,t)| < e/V(f[Xc]. We define 

F{uj) = r^iJ,t)X{dt) =J2dxn{xi,X)fix{{t : iu,t) G d}), 
•'^ i 

where fix is the Lcbcsguc-Stieltjes measure on [0, T] corresponding to 
the function A. Fubini's theorem implies that the function uj IJ-xiit '■ 
{io, t) G Ci}) is C7c-measurable. The inequality (12) is obviously satisfied. 

□ 

Proof of Theorem 2. We take a sequence of partitions = {0 = to < 
< ■ ■ ■ < = T}, and consider the functions 

n-l 

x^{t) = j2xm^,.,t^^)it) + xiT). 

i=0 



Noneqworkthm.tex; 31/03/2008; 15:39; p. 8 



9 



Clearly, A„ ^ A on [0,r], and V^[X^] < V^[\], where is the 
variation on [0,T]. We set A" = A(t"), and define the functions 

^n{^) = e-/^E:=i(WM*?).^?)->^K*?).^"-i)). 

By Theorem 1, 



/ 



(Pn{uj)L,xn{duj) - 



I Til' \— / — /\ V / ry 

/n ^A{o) 
for all n. We denote (p{uj) = e~^^^^'^\ and prove that 

lim / {(pn{(^) — ip(u!))lL\n(du!) = 0, (13) 

lim / {ifniuj) - (^(w))LA((iw) = 0, (14) 
"^<» 

lim lim / (<z)„(a;) — (/?(a;))L\m(cL!) = 0. (15) 

n— >oom— +00 Jq 

For this, we first replace the functions cpn in (13)-(15) with more suit- 
able functions ifn such that (99^(0;) — (^^(a;)) ^ as n — >■ 00, uniformly 
on n. We have: 

f:{nHt2), A?) - n{u;{t2), K-i)) - f d^n{u{t), \^{t)) dx^it) 

n n 

= ^ W(a;(tr), Ar)(A? - Ar_i) - ^ W(c.(C), A?)(Ar - A^^J, 

i=l 1=1 

(16) 

where in the first term in (16) we applied the mean value theorem 
to each summand, and chose A" G [Af , A"_,_]^]. Since by assumption 
dxH{x, ■ ) is equicontinuous, the absolute value of the difference in 
(16) does not exceed eVg^[A] where e is chosen so that \dxTC{x,Xi) — 
d\H{x, \2)\ < e whenever |Ai — A2I < 5, and 6 is chosen by the equicon- 
tinuity argument. The relation (16) shows that if we prove (13)-(15) 
with <p„{uj) = e-f^Iodxn{^{t),xr^{t))dX-{t) substituted for cp„, then we 
prove (13)-(15). We define the functions: 

a.-(t,a.) =a;(0) + Er=i'^(ir)%_„*r]W; (17) 

A«(t)=A(0)+E?=lA(i?)I(in_^,tn](i). 
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With the help of these functions, the second term in (16) can be 
represented as 

r dxn{u;{t),X'\t))dy{t) = r dxH{u''{t,uj),X'{t))dX^{t) 
•^0 -^Or (18) 

Jo 

By Assumption 7, we choose an e > so that |(?A'H(w"(t, w), A"(i)) — 
dx7i{oj'^{t,ui),X(t))\ < £ whenever sup^g^o,!] I^^C*) ~ < ^i ^"^^ ^ is 
chosen by the equicontinuity argument. This means that the relations 
(13), (14), and (15) are equivalent to 

hm / {(pnioj) — u)iuj))'h\n{dijj) = 0, 

hm / {0n{to) - ip{u;))h),{du;) = 0, (19) 

lim hm / ((5„(a;) — ip(u))h\m (du) = 0, 

n— >oo m^oo 

where (^„(u;) = e-^^/cr 9AW{<i"{t,<^),A(t))dA(t) ^ ^"(t,a;) is given by (17). 
We show that the relations (19) hold. Since dxH{ ■ , ■) is bounded by 
Assumption 6, A is a function of bounded variation on [0, T], and the 
exponent is Lipschitz on bounded domains, for all m we obtain the 
estimate 

((^n('^) - ¥'(a;))LAm((ia;) (20) 

<Kl [ ILxAduj) r \\\{dt)\dxn{u''{t,Lo),\{t)) - dxH{Lo{t),\{t))\ 
Jn Jo 

= Kl r\X\idt) I l.xAdu)\dxn[Cj'\t,uj),\{t)) - dxH{u:{t),\{t)% 



where is the Lipschitz constant for the exponent, |A|((it) is the 
Lebesgue-Stieltjes measure corresponding to the total variation func- 
tion |A|(t), and Fubini's theorem has been applied to pass to the last 
integral. The same estimate holds for L;^. Namely, 



{0n{^) - ^{^))^\{duj) 

n 

T 



(21) 

Kl r \X\{dt) [ hx{dw)\dxn{u;"{t,u;),X{t))-dxn{Lo{t),Xm- 
Jo Jn 

We would like to show that 

lim / lLxr^(duj)\dxn{u''\t,uj),X{t)) -dxn{oj{t),X{t))\=0, (22) 
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X 



lim / l.x{di^)\dxH{Cj'\t,u),X{t)) -dxH{u{t),X{t))\=Q, (23) 
lim lim / LAm(cL;)|5AH(a;"(t,a;),A(t)) -aAH(a;(i),Am)|=0. (24) 

n— >oo m— »oo Jj-^ 

Let t G (i"_i, i"], and let 8n = t'^ - t. Then a;"(t, w) = cc;(t + (5„). Let us 
show (22). We have: 

hxnidco) \dxH{u{t + 5n),X{t)) - dxn{uj{t),X{t))\ 

qXo{xo)dxo / dxpxn{to,xo,t,x) 

JX 

I dxip{t,x,t + 6n,Xi,X{t)) \dx'H{x,X{t)) - dxHixi, X{t))\. 

JX 

Note that 

Fn{x) = / dxip{t,x,t + dn,Xi,X{t)) dxH{xi,X(t)) -dxT-i{x,X{t)) 

JX 

converges to zero uniformly in x running over compacts in X which 
follows from Assumption 3. Indeed, 

|F„(x)|2^ / dxip{t,x,t + dn,Xi,X{t)) dxH{xi,X{t))-dxn{x,X{t))^. 

JX 

The right hand side of this inequality converges to zero uniformly in 
x K C X, where K is an arbitrary compact. This easily follows 
from Assumption 3 after we separate in the right hand side the terms 
depending on x and on a;,. Define the measures: 

/^A»(-4) = Jxq\o{xo)dxo J^pxr^{to,xo,t,x)dx, 
/^a(^) = Jx Q\)(xo) dxo J^px{to,xo,t,x) dx. 

By Assumption 4, L^" converges weakly to relative to the family of 
bounded continuous cylinder function, whenever A" ^ X. This implies 
that nx^ nx weakly (relative to the family of bounded continuous 
functions), as A" ^ X. By Prokhorov's theorem, the family {/J-x^, fJ'x} of 
probability measures on X is tight. We fix an arbitrary e > and find 
a compact such that fix"iX\Ki;) < e for all n. Since F„ ^ on 
Kf,, and all Fn are bounded on X by a constant, say M, not depending 
on n, we can find an A'' G N such that |F„(x)| < £ on ATg for all n > N. 
We obtain the estimate: 

/ lxxnidx)Fnix) ^ sup \Fn{x)\ + Mnxr^{X\Ke) ^e + Me. 

JX xeKe 

This proves (22). The relation (23) follows from the right continuity 
of Lo, Assumption 8, and Lebesgue's theorem. The relation (24) follows 
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from Assumption 4, and from (23). Application of Lebesgue's theorem 
to the integral (20) (taken over [0,T] with respect to |A|((it)) imphes 
that as n ^ oo, the integral (20) converges to zero. This implies (19), 
and thus (13), (14), and (15) are proved. By (13), we obtain: 

= lim / ((pnLixn — ipl^\n) = lim / ipnLxn — lim / (fhxn 

n— >oo J n— >oo J n— ►oo J 

lim / ipLixn. 
n— »oo J 



Z\(T) 



On the other hand, (14) and (15) imply: 

= lim lim / {U)\nl^xm — iphxm) 

= lim lim / (pnLxm — lim / iphxm 

n— >oo m^oo J m— »oo J 

= lim / ipn^x — lim / ipl^x"" = / ~ lim / V'La"'- 
n— >oo J m— >oo J J m— »oo J 

Comparing the last two relations gives: 

/ ^{u;)hxi<ko) = 

^A{0) 

The theorem is proved for the case A G C^[0,r]. □ 



3.3. JARZYNSKI'S identity for a G F[0,r] AND ITS COROLLARIES 

Corollary 1 (Corollary of Theorem 2). Let the assumptions of The- 
orem 2 be fulfilled, and let f : X — > M 6e bounded and continuous. 
Then, 

EL,[(/o7rT)e-^^^] =E,(^)[/]ELje-^^^], (25) 

where n : X^'^l ^ X, u 1-^ Lo(t) is the evaluation mapping, ^x{T) is 
the expectation relative to the measure — e^^^^^'^^-^^^dx. 

Proof. Assuming that A G /^stci)[0; ^]) '^^ repeat the argument of (6) 
and (10), while using the relation (9). Specifically, we obtain: 

EL,[(/o7rT)e-/5^^] = / e-^«(--^(*"))/(rr„)dx„ 

^A(to) J. 

f p-/3H(x„_i,A(t„_i)) 

■ P{tn—li ^n—li tfii ^ni -^(^n— l)) dx^—l . . . 



g-/3W(a:„,A(t„_i)) 



X 
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l3H{xo,Xito)) 



A(to) 

Now let A G C^[0, T], and let the constant Mf be such that 
sup^^x 1/(^)1 < -^/- -^s in the proof of Theorem 2, we find a sequence 
of step functions A'* converging to A uniformly on [0, T]. The identity 
(25) holds for each A". The argument of passing to the limit as n — > oo 
is similar to the argument used in the proof of Theorem 2. We define 
the function F : O — > M, F{lo) = f{uj{T)), and repeat the argument of 
Theorem 2 until the inequality (21) with the following replacements: 
(p Fif, ifn ^ Fipn, (pn ^ F(pn, 0n ^ Fipn. Thc right hand sides 
of (20) and (21) will be the same as in the proof of Theorem 2 but 
the Lipschitz constant Kl will be replaced with K^Mf. The part of 
the proof of Theorem 2 following after the inequality (21) remains 
unchanged until the last two arguments of passing to the limit. Those 
arguments now will be: 

= lim / {F(pnL,xn — F(pL,xn ) = lim / FiprX'X^ ~ 1™ / F(pL,xn 

n— >oo J n^oo J n— >oo J 

= ^x{T)[f] ELje-^^^] - Jim j FypLAn. 
On the other hand, 

= lim lim I (Fipxnl^x'^ ~ Fiphx"^) 
n— >oom— »oo J 

= lim lim / F'x>n'L,xm — lim / Fiphx^ 

ra— >oom— >oo J m— >oo J 

= j F<^La -^im^ j Fiphx^- (26) 
This implies the identity (25). □ 



Corollary 2 (Corollary of Theorem 2). Let the assumptions of Corol- 
lary 1 except Assumption 4 of Theorem 2 be fulfilled with respect to the 
function A G C^[0,T]. Let a E A, and let A" = AI[o_t) + ^^{t}- Let us 
assume that 

Ify e Cstcp[0,T], and A" ^ A" on [0,T], then La- weakly 

relative to the family of bounded continuous cylinder functions; 
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Then, 

EL,a[(/o7rT)e-^^^] =E,(^)[/] (27) 

Proof. If A« G £step[0,r], i.e. A" = E^j/ >^^h^,ti+l) + « V}' t^ien the 
representation (4) shows that El^^o [(/ o ttt) e"/^^^] is the same for the 
functions A" with different a. Hence, 

lEL.a [(/ o ttt) e-^ = E,(r) [/] ^ . 

Note that in the relations (20) and (21) we integrate with respect to 
|A|(dt) where |A| is continuous. Hence, by Lebesgue's theorem, the 
investigation of convergence at the point T in the integrals with re- 
spect to h\^{dio) and hx{dio) is not necessary. Finally we note that by 
assumption, and by the relations (21) and (23), with A" substituted for 
A, for the first term in the identity (26) we obtain: 

lim lim / Fip^hx^ = / F(p'L\a. 

n— ►oom— »oo J J 

This proves (27). □ 

Let [s,t] C [0, T] be an arbitrary subinterval, Xl'*'*! be the space of 
all functions on [s, t]. We introduce further notation: Let C^[s, t] be the 
space of continuous functions of bounded variation on [s,t], jCstep[s,t] 
be the space of right continuous step functions on [s, t] taking a finite 
number of values. Below, on Xl^'*! we define the distributions L^.^ ^ 
and I^\;s,t, A G V[0,T], a; G X. Let s < ri < ■ ■ ■ < Tk-i < t he a. 
partition of {s,t\, f : X^~^^ ^ M be bounded and measurable. The 
finite dimensional distributions of L^.^ ^ are defined by the formula: 

f{uj{s),uj{Ti), uj{t)) L^;,,j (dw) = 
j dxi px{s, X, ri, xi)) ■■■ j dxkPx{Tk-i,Xk-i,t, Xk) f{x, xi,..., Xk). 

X X 

We extend ^x-st o'c(^^^'*') by Kolmogorov's extension theorem, and 
notice that under assumptions that we made, this measure it concen- 
trated on the right continuous paths without discontinuities of the 
second kind. Let / be L^.^ ^-integrable. We define 

^X;sAf] =^X{s)Hs,M- (28) 
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Theorem 3 (Jarzynski's identity, case: A G V[0,T]). Let A G V[0,T], 
and let {ti < . . . < tk = T} be the set of its discontinuity points. 
Further let the transition density function p\ of the Markov process 
satisfy Assumptions 1, 2, and 3, and the probability distribution hx be 
given by (2). Let Assumptions 5, 6, 7, and 8 of Theorem 2 be fulfilled. 
Additionally, we assume that 

8. If e Cstep[0,T], and ^ \ on [0,r], as m ^ oo, then for 
each i, ^ i < n, l^\'n-ti,ti+i ~^ ^X;ti,ti+i weakly relative to the 
family of bounded continuous cylinder functions on 

9. For all bounded and continuous functions f : — > M, the 
function X — > M, x i— > ^A-tj u+i if] continuous. 

Then, the function e~^^^ is hx-integrable, and 

ElJc-'^^-] = e-^^^. (29) 

The continuity of functions — > M is understood with respect 
to the topology of pointwise convergence. 

Proof of Theorem 3, case A G y[0,T]. Let A = Ac + Agtep, where Ac G 
C^[0,r], and Astep = T.'lZ^ h\u,U+.) + e >Cstep[0,r]. On each in- 
terval we define the function uii = and identify each 

u G Xt'^'-^l with the sequence (wi, . . . , w^). Since is a Markov process, 
we obtain: 

g-/3WA(a;i,...,a;„) 



n „f. 

Wxicvi,...,un) =T. I dxniiv^{t),Xit))dX,{t) 



where 

n „f. 
n 

+ Y^iuiuiiti), \{ti)) - n{ui{ti),\{ti - 0))). 
1=1 

On each interval we define continuous functions \{i){t) = A(t), 

t G \ti-i,ti), and A(j)(ii) = \{ti — 0), 1 ^ i ^ n, and introduce a 
notation: 

Wx^^{ui)= r dxn{ui{t),x{t))dx^i^{t). 
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We obtain: 

-m. 



■^(i)('^l)e-/3(W(c^i(tl),A(ti))-W(a;i(ti),A(ti-0)))La;o ,g \ 



X['l'*2] 



/ A;t„_i,t„ V ny- 

Xl*n-l>*nl 

Let denote the expectation relative to the measure q^{x)dx. We can 
rewrite the above relation in terms of expectations: 



(30) 



Let F : ^ M be bounded and continuous. We obtain: 



E 



'A(ti-O) 



g-/3(H(.,A(ti))-H(.,A(ti-0)) 



-^A(ti-O) 



E 



■^A(ti-O) 



(31) 



Note that in (30) F is always a function of the form 

^ M, u; ^ e-"^'^(^+-^)("V(^(t.+i)), 



(32) 



where / : X ^ M is bounded and continuous by Assumptions 8 and 
9. Taking the expectati 
Corollary 2 we obtain: 



9. Taking the expectation El;^.^ of the function (32), and applying 



E- 



[e-^^^..i)('^)/(a;(t,^,))] 



Z 



Applying the identities (31) and (33) to (30), we obtain: 



EA(t,+,-o)[/]- (33) 



^Xitn-l) ^A(tn-O) 



-A(T) 
^A(O) 



• (34) 



□ 
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Corollary 3 (Corollary of Theorem 3). Let the assumptions of The- 
orem 3 he fulfilled, and let f : X — > R 6e bounded and continuous. 
Then, 

ELj(/o7rr)e-'^'^^] =E,(^)[/] ELje"^^^]. (35) 

Proof. Wc repeat the argument that we used in the proof of Theorem 
3 to obtain (30) in connection to the expression El;^[(/ o ttt) e^'^''^^]. 
Instead of the very last expectation in (30), we obtain 

\-„_i(t„_i)[e"^^^")^""^ y(^^(^^)^g-/3(H(a;„(tn),A(U))-W(c.„(t„),A(t„-0)))]_ 

Applying Corollary 2, we obtain: 

^ ^A(t„-0) ■ Mtn))-n( ■ Mtn-O))) J] 

_^A(t„-0) ^A(i„) (36) 



^A(t„_i) ^A(t„-0) 



z z 

Instead of the multiplier ^ in the relation (34) we obtain 



2^7^ 7 zVir^ ^A(t„)[/] according to (36). This gives (35). □ 



4. Bochkov— Kuzovlev's identity 

Here we give a rigorous mathematical proof of the identity announced 
in [4]-[7] where the authors used a different definition of work to that 
in the papers [l]-[2]. The difference between the two definitions of work 
was analyzed in the paper [3] , and it was found that 

W{u) - Woiio) = n{u{T), A(T)) - n{co{T), A(0)), 

where Wq is the work in the Bochkov-Kuzovlev sense. We use this 
equality as the definition of the new work Wq, and will prove the 
following theorem: 

Theorem 4. Let the assumptions of Theorem 3 he fulfilled. Then, 

EL,[e-^'^°] = l. (37) 
Moreover, the identities (37) and (29) are equivalent. 
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Proof. Applying Corollary 3 and the identity (29) we obtain: 

This relation also shows that the identities (37) and (29) are equivalent. 

□ 



Acknowledgements 

The author thanks Professor C. Dellago for attracting her attention 
to Jarzynski's identity, and Professor O. G. Smolyanov for useful dis- 
cussions. This work was supported by the research grant of the Erwin 
Schrodinger Institute for mathematical physics, by the Austrian Sci- 
ence Fund (FWF) under START-prize-grant Y328, and by the research 
grant of the Portuguese Foundation for Science and Technology (FCT). 



References 

1. C. Jarzynski, Nonequilibrium equality for free energy differences, Phys Rev 
Lett 78, 2690, pp. 2690-2693, 1997. 

2. C. Jarzynski, Equilibrium froo-oncrgy differences from nonequilibrium mear 
surements: A master-equation approach, Phys. Rev. E 56, 5018-5035, 1997. 

3. C. Jarzynski, Comparison of far-from-equilibrium work relations, Comptes 
Rcndus Physique, 8, 495 506, 2007. 

4. G. N. Bochkov, Yu. E. Kuzovlev, General theory of thermal fluctuations in 
nonUnear systems, Zh. Eksp. Teor. Fiz. 72, 238-243, 1977; [JETP 45, 125-130, 
1977]. 

5. G. N. Bochkov, Yu. E. Kuzovlev, Fluctuation-dissipation relations for nonequi- 
librium processes in open systems, Zh. Eksp. Teor. Fiz. 76, 1071-1079, 1979; 
[JETP 49, 543-551, 1979]. 

6. G. N. Bochkov, Yu. E. Kuzovlev, Nonlinear fluctuation-dissipation relations 
and stochastic models in nonequilibrium thermodynamics, I. Generalized 
fluctuation-dissipation theorem, Physica 106A, 443-479, 1981. 

7. G. N. Bochkov, Yu. E. Kuzovlev, Nonlinear fluctuation-dissipation relations 
and stochastic models in nonequilibrium thermodynamics, II. Kinetic potential 
and variational principles for nonlinear irreversible processes, Physica 106A, 
480, 1981. 

8. G. E. Crooks, Path-ensemble averages in systems driven far from equilibrium. 
Physical review E 61, pp. 2361-2366, 2000. 

9. G. E. Crooks, Nonequilibrium Measurements of Free Energy. Differences for 
Microscopically Reversible Markovian Systems. Journal of Statistical Physics, 
90, Nos 5/6, pp. 1481-1487, 1998. 



Noneqworkthm.tex; 31/03/2008; 15:39; p. 18 



19 



10. G. E. Crooks, Entropy production fluctuation theorem and the nonequilibrium 
work relation for free energy differences, Physical Review E 60, pp. 2721-2726, 
1999. 

11. C. Dellago, E. Scholl-Paschinger, A proof of Jarzynski's non-equilibrium work 
theorem for dynamical systems that conserve the canonical distribution, J. 
Chem. Phys. 125, 2006. 

12. C. DcUago, P. L. Gcisslcr, W. Lcchncr, H. Obcrhofcr, Equilibrium free ener- 
gies from fast-switching simulations with larg time steps. Journal of Chemical 
Physics 124, 044113, 2006. 

13. E. B. Dynkin, Markov processes. Academic press, 1965. 

14. G. Hummer, A. Szabo, Free energy reconstruction from nonequilibrium single- 
molecule pulling experiments, Proc. Nat. Acad. Sci. 98, 3658, 2001. 

15. G. Hummer, A. Szabo, Free energy surfaces from single-molecule force 
spectroscopy. Acc Chem Res 38, 504-513, 2005. 

16. J. Liphardt ct al.. Equilibrium information from nonequilibrium measurements 
in an experimental test of Jarzynski's equality. Science 296, 1832, 2002. 

17. N. I. Portenko, A. V. Skorokhod, B. M. Shurenkov, Markov processes, 
Itogi nauki i techniki. Sovremennye problemy matematiki. Fundamentalnye 
napravleniya. 46, pp. 5-248, 1989 (in Russian). 



Noneqworkthm.tex; 31/03/2008; 15:39; p. 19 



Noneqworkthm.tex; 31/03/2008; 15:39; p. 20 



